subroutine comput_flux2
use var
implicit none

    integer m,j,i
    
        do m=1,4
            call flux2_m(m)
            call flux2_p(m)
            do i=1,gridsx
                do j=0,gridsy
                    flux2(i,j,m)=fp2(i,j,m)+fm2(i,j,m)
                enddo
            enddo            
        enddo
        
end subroutine